TO BE PUBLISHED THE GALAXY EVOLUTION EXPLORER (GALEX) ASTROPHYSICAL JOURNAL LETTERS SPECIAL ISSUE 
Preprint typeset using LJTgX style emulateapj v. 4/9/03 



GALEX OBSERVATIONS OF THE ULTRAVIOLET HALOS OF NGC 253 AND M82 

Charles G. Hoopes 1 , Timothy M. Heckman 1 , David K. Strickland', Mark Seibert 2 , Barry F. Madore 34 , R. Michael 
Rich 5 , Luciana Bianchi 6 , Armando Gil de Paz 3 4 , Denis Burgarella 7 , David A. Thilker 1 Peter G. Friedman 2 , Tom A. 
Barlow 2 , Yong-Ik Byun 8 , Jose Donas 7 , Karl Forster 2 , Patrick N. Jelinsky 9 , Young- Wook Lee 8 , Roger F. Malina 7 , D. 
Christopher Martin 2 , Bruno Milliard 7 , Patrick F. Morrissey 2 , Susan G. Neff 10 , David Schiminovich 2 , Oswald H. W. 
Siegmund 9 , Todd Small 2 , Alex S. Szalay 1 , Barry Y. Welsh 9 , and Ted K. Wyder 2 

Draft version February 2, 2008 



o 
o 

(N 
> 

o 



> 

OS 

o 



ABSTRACT 

We present Galaxy Evolution Explorer (GALEX) images of the prototypical edge-on starburst galaxies M82 
and NGC 253. Our initial analysis is restricted to the complex of ultraviolet (UV) filaments in the starburst- 
driven outflows in the galaxy halos. The UV luminosities in the halo are too high to be provided by continuum 
and line emission from shock-heated or photoionized gas except perhaps in the brightest filaments in M82, 
suggesting that most of the UV light is the stellar continuum of the starburst scattered into our line of sight by 
dust in the outflow. This interpretation agrees with previous results from optical imaging polarimetry in M82. 
The observed luminosity of the halo UV light is < 0.1% of the bolometric luminosity of the starburst. The 
morphology of the UV filaments in both galaxies shows a high degree of spatial correlation with Ha and X-ray 
emission. This indicates that these outflows contain cold gas and dust, some of which may be vented into the 
intergalactic medium (IGM). UV light is seen in the " Ha cap" 1 1 kpc North of M82. If this cap is a result of 
the wind fluid running into a pre-existing gas cloud, the gas cloud contains dust and is not primordial in nature 
but was probably stripped from M82 or M8 1 . If starburst winds efficiently expel dust into the IGM, this could 
have significant consequences for the observation of cosmologically distant objects. 

Subject headings: galaxies: halos — galaxies: starburst — galaxies: individual(M82, NGC 253) — ISM: jets 
and outflows — ultraviolet: galaxies 
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1. introduction 

Many local starburst galaxies have galactic-scale outflows 
of metal-enriched gas, called starburst superwinds, which are 
driven by the stellar winds and supernovae of numerous mas- 
sive stars (e.g., Heckman, Armus, & Miley 1990). These out- 
flows contain a hot (10 7 K), met al-enriched w ind fluid as wel l 
as entrained cooler gas and dust ( Strickland & Stevens 2000). 
It is likely that the hot gas can escape from the potential 
well of the parent galaxy, enrich ing the i ntergalactic medium 
(IGM) with metals and energy (Heck man et alJ l2000). It is 
not clear whether the colder gas can escape, and this is an 
important question since it would mean that superwinds also 
enrich the IGM with du st, which coul d affect observations 
of high-redshift objects |Aguirrdll999t lAguirre et al.ll2001l 
lAlton. Davies. & Bianchilll999t iHeckman et alJl2000l . This 
question is even more crucial since similar outfl ows are now 
known to be common in high-redshift starbursts (tPettini et al.l 
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l27MISharTev et al.l2003l) . 

Several lines of evidence suggest that su perwind s con- 
tain d ust iHeckman et al.l l2000 : Shonbell & Bl and-Hawfhornl 
1998). Optical imaging polarimetry shows light scattered by 
dust in the halos of starburst galaxies, including M82 (e.g 
Scarrott, Eaton, & Axon 1991; Alton et al. 1994). Far-IR 
and sub-mm imaging reveal thermal emission from extrapla- 
nar dust in several edge-on starburst galaxies (Alton, Davies, 
& Bianchi 1999). Finally, the strong correlation between 
the strength of the blueshifted interstell ar Na D line a nd the 
line-of-sight reddening in superwinds (Heck man et al.1 2000). 
strongly suggests the dust is actually outflowing. What re- 
mains unclear is the physical relationship between the cool 
dust-bearing gas and the warm and hot gas probed by optical 
lines and X-rays respectively. Comparison of sensitive, high 
resolution images of the dusty material with Ha and X-ray 
emission images would shed light on this relationship. Dust 
is highly reflective in the ultraviolet (UV) (Draine 20031), so 
imaging of starburst superwinds in the UV can trace the lo- 
cation of dust, if one can account for UV emission by pho- 
toionized or shock-heated gas. Indeed, Ultraviolet Imaging 
Telescope (UIT) near-UV data for M82 show evidence of 
UV light in the halo corresponding to known Ha features 
( Marc um et al I200T1) . Here we present Galaxy Evolution Ex- 
plorer (GALEX) ultraviolet images of two prototypical star- 
burst superwind galaxies: NGC 253 and M82. The images 
reveal prominent UV light in the superwind region. Our goal 
is to understand the origin of this light. 

2. OBSERVATIONS 

NGC 253 was observed by GALEX on 2003 October 13 
for 3289 seconds. M82 was observed by GALEX on 2003 
December 8 for 3083 seconds. The GALEX data include 
far-ultraviolet (FUV; X eff = 1528 A, AA = 268 A) and near- 
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FIG. 1. — NGC 253 in UV and Ha. The left panel shows a two-color image, with GALEX NUV in red and FUV in blue. The right panel shows the 
continuum-subtracted Ha image. The images are 30' (22.7 kpc) on each side, with North up and East on the left, and are aligned with each other. The intensity 
scales in both panels are logarithmically scaled and stretched to emphasize the faint, diffuse emission, so the bright disk of the galaxy is saturated. 



ultraviolet (NUV; X eff = 221 \ A, AA = 732 A) images with a 
circular field of view with radius ~ 38'. The spatial resolution 
is ~ 5 ". Details of the GALE X instrumen t and dat a character- 
istics can be found in Martin et al. ( 2004) and Morrissev et al. 
(2004). 

We also use previously obtained Ha data. The 
Ha image of NGC 253 is desc ribed in detail in 
iHoopes. Walterbos. & Greenawaltl ( fi99r3) . The Ha image of 
M82 is part of a mosaic of the M81-M82 system taken with 
the Burrell-Schmidt telescope at KPNO, and is described in 
lGreenawaltet "aTI Jl998h . 

3. ANALYSIS 

3.1. Ultraviolet Morphology 

Figure ^ compares the two-color GALEX image of 
NGC 25 3 with the Ha image. Extended Ha emission was 
noted by IStrickland etafl |2002). Diffuse emission extend- 
ing several kpc into the halo on both sides of the disk, north- 
east of the galaxy center with the brightest an d more ext ended 
emissi on toward the east end of the disk. Stricklan d et alJ 
(2002) found that the X -ray emission matched the Ha emis- 
sion in morphology. These features are also visible in the 
GALEX images. 

Figure|2]shows the UV and Ha images of M82 (in the same 
manner as Figure^. The M82 images show a bright, com- 
plex network of filaments, very different in appearance from 
NGC 253. The morphology in the UV and Ha images is 
strikingly similar. Prominent Ha filaments are seen perpen- 
dicular to the disk on both sides, surrounded by a lower sur- 
face brightness component of diffuse light (see also Ohyama 
et al. 2002). The filaments are also visible in the UV, but 
there is less contrast between the filaments and the diffuse 
UV light. The GALEX image s are much more se nsitive than 
the earlier UIT NUV image dMarcum et alJl200ll) . and show 
that the UV- Ha correlation in morphology extends to very 
faint Ha filaments, and exists in the FUV as w ell (which 
was not detected by UIT). IStrickland et alJ J2004I) noted that 
the X-ray and Ha morphology are similar on all scales, and 



this is also true for the UV light. The Ha and X-ray "cap" 
JDevine & BallvlH999HLehnert et aljfl999h 1 1 kpc above the 
North side of the disk is visible in both GALEX bands. We 
will address the UV properties of these and other nearby star- 
bursts, including possible reasons for the striking differences 
between M82 and NGC 253, in a forthcoming paper (Hoopes 
et al., in preparation). 

3.2. Luminosities and Flux Ratios 

Table^compares the UV and Ha luminosities of the halo 
with the total and bolometric luminosities. The measurements 
have been corrected for Galactic foreground extinction using 
E(B-V) = 0.019 for NGC 253 and E(B-V) = 0.493 for M82 
JSchlegel. Finkbeiner. & Davislll998t) . A correction factor of 
0.59 has been applied to remove [N II] from the Ha flux. The 
extraplanar UV light in both cases is less than 0.1% of the 
bolometric luminosity of the starburst. The observed halo lu- 
minosity is 7% (10%) of the total observed NUV(FUV) lumi- 
nosity of NGC 253, and for M82 it is 43% (65%). The Ha 
luminosity of the halo is 4% of the total Ha luminosity for 
NGC 253 and 2 1 % for M82. 

Figure [5] shows flux ratios measured in square regions 
30" on each side. The GALEX monochromatic fluxes were 
multiplied by the effective filter bandpass to give units of 
ergs cm" 2 s . Figure [3] also shows model predictions for 
the continuum (Balmer, Bremsstrahlung, and two-photon) 
and line emission of shock-heated and photoionized gas 
( Dopit a & Sutherland! 1 9961 iFerlandl 1996ft . The shock mod- 
els span shock velocities from 100 to 900 km s" 1 , and include 
both shock and precursor emission. The photoionization mod- 
els are of spherically symmetric clouds ionized by a central 
source (the UV continuum of the ionizing source is not in- 
cluded in the model predictions), and span stellar tempera- 
tures ranging from 30000 K to 50000 K and electron densities 
from 0.1 to 10 cm" 3 . Solar abundances were assumed in both 
cases. 

Most of the regions have too much UV light (relative to 
Ha) to be explained by nebular emission alone. The observed 
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FIG. 2. — M82 in UV and Ha. The left panel shows a two-color image, with GALEX NUV in red and FUV in blue. The right panel shows the continuum- 
subtracted Ha image. The images are 21' (22.0 kpc) on each side, with North up and East on the left, and are aligned with each other. The intensity scales in 
both panels are logarithmically scaled and stretched to emphasize the faint, diffuse emission, so the bright disk of the galaxy is saturated. 



Table 1. Measured Properties" 



Galaxy 


Distance 


Lua Halo 


La a Total 


Lnuv Halo 


Lnuv Total 


Lfuy Halo 


Lfuv Total 


L Bo l b 






(Mpc) 


(erg s" 1 ) 


(erg s~') 


(erg s -1 ) 


(erg s" 1 ) 


(erg s~' ) 


(erg s~') 


(erg s~') 




NGC 253 


2.6 


1.5 x 10 39 


3.8 X 10 40 


3.1 X 10 40 


4.5 X 10 41 


2.1 X 10 40 


2.2 X 10 41 


7.8 X 10 43 


-1.5 


M82 


3.6 


1.3 x 10 40 


6.1 x 10 40 


1.5 X 10 41 


3.5 X 10 41 


7.1 X 10 40 


1.1 X 10 41 


2.0 X 10 44 


-0.6 



a The measured luminosities were corrected for foreground Galactic extinction. Calibration uncertainties are ~ 10 % in the UV bands iMorrissey et al. 2004) 
and are of similar magnitude in Ha. 
b The bolometric luminosities were taken from lStrickland et all OoTjOl) (NGC 253) and lMcLeod et alJ<T^3l) (M82). 
c The spectral slope, defined via F\ oc A' 3 , was estimated from the FUV/NUV flux ratio following Kong et al. 1 2004). 



FUV/ Ha ratios in some of the brighter regions in the M82 
halo are consistent with a significant contribution from shock 
ionization, but these values have not been corrected for ex- 
tinction intrinsic to the wind. The optical spectrum of the 
M82 wind indicates a reddening of E(B-V) > 0.21 mag 
jHeckman et al.lH990l) . which would increase the observed 
FUV/Ha ratios by a factor of > 3.9 (> 2.6 for NUV/Ha). 
This implies that the wind is substantially brighter in FUV 
and NUV than would be possible for photoionized or shock 
heated gas. The absence of O VI emission seen in Far Ultra- 
violet Spectroscopic Explorer data limit s shock speeds in th e 
bright M82 filaments to v s < 160 km s" 1 (Hoopes et al. 2003), 
much slower than the wind velocity (v > 10 3 km s" 1 , Strick- 
land & Stevens 2000). Taken together, these facts imply that 
another source is required to explain the excess extraplanar 
UV light. The diffuse morphology argues against star forma- 
tion in the wind as the source. The most likely remaining 
mechanism is scattering of stellar continuum from the star- 
burst by dust in the halo. 

4. DISCUSSION 

The spectral slopes (3 in the halo implied by the observed 
FUV/NUV ratio are listed in Table ^ The values are redder 
than an unreddened starburst (-2.5 < (3 < -2.0), and in fact 
agree well with observed (/. e., reddened) values of local star- 
bursts (-2.0 < (3 < -0.6; Meurer, Heckman, & Calzetti 1999). 



While it is clear that dust in a starburst environment may have 
pro perties that differ from the standard Galactic dust m od- 
els dGordon . Calzetti, & Wittll 19971 iPopescu et aTll2000l) . the 
FUV/NUV ratio i s in general agreement with the dust scat- 
tering models of (Draine 2003) in which the dust albedo is 
greater in the NUV than in the FUV. 

In a sample of local star forming galaxies. lBuat et aflj2002l) 
found that the observed Ha flux was on average 2.5-5.0% of 
the observed UV flux near 2000 A (similar to the GALEX 
NUV band). The observed Ha flux of the NGC 253 halo is 
5% of the NUV flux, and for M82 the corresponding value 
is 12%. This may indicate that not all of the Ha emission 
in the M82 halo is scattered light, although polarization mea- 
surements indicate the presence of some scattered light in the 
Ha filaments and the more extended diffuse Ha component 
of M82 (e. g., Scarrott et al. 1991). Scattered light could be a 
substantia l component of the fainter NGC 253 halo. 

iLehnert et alJ ( 1 999 ) suggested that the Ha cap near M82 is 
the result of a collision between the hot wind fluid and a pre- 
existing neutral cloud. If this scenario is correct, our results 
imply that the cloud is not primordial since it contains dust. 
The cloud may have been pushed out of M82 by the starburst 
wind or stripped from either M82 or M8 1 by the tidal interac- 
tion between the two galaxies. 

Our results establish for the first time a close morpholog- 
ical correspondence between the dust and the hotter phases 
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FIG. 3. — Comparison of measured ratios with model predictions. The points are the measured values, and have been corrected for Galactic foreground 
extinction using the extinction law of Cardelli, Clavton, & Mathis 1 1989). The hatched regions indicate the range of model predictions. The models are described 
in the text. We have not attempted to correct for internal extinction. Reddening vectors for E(B-V) = 0.21 (measured in the M82 wind by Heckman et al. 1990) 
are shown in each panel, assuming the Calzetti 1 2001) starburst extinction law. The M82 points are separated into bright and faint based on their Ha surface 
brightness, with the division occurring at I(Ha) = 4 X 10~ 17 ergs s _1 cirT 2 arcsec -2 . The two circled triangles are in the M82 Ha cap. Representative error bars 
are shown in each panel. 



of the winds probed in Ha and X-ray emission. We have 
direct evidence that the hotter gas is outflowing (e.g., Strick- 
land & Stevens 2000), so the new UV images provid e further 
evidence for outflowing dust (Heckma n et alJl2000l) . This is 
consistent with the idea that the cool dusty material is ambient 
interstellar gas in the disk or inner halo that has been has been 
entrained and accelerated by the hot outflowing gas generated 
in the starburst. 

If this dust is ejected into the IGM, there could be important 
implications for cosmological observations. While the dust 
density is small, over cosmological dista nces t he resulting ex- 
tinction could be sig nificant (Aguirre 1999; Heckm an et al.l 
2000). lAlton. Bianchi. & Daviesl (120011) point out that inter- 
galactic dust may affect the determination of the evolution of 
the cosmic star formation rate, for example. More work is 
needed to understand the effects of intergalactic dust. 
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